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We derive analytical expressions for the conductivity of bilayer graphene (BLG) using the Boltz- 
mann approach within the the Born approximation for a model of Gaussian disorders describing 
both short- and long-range impurity scattering. The range of validity of the Born approximation 
is established by comparing the analytical results to exact tight-binding numerical calculations. A 
comparison of the obtained density dependencies of the conductivity with experimental data shows 
that the BLG samples investigated experimentally so far are in the quantum scattering regime where 
the Fermi wavelength exceeds the effective impurity range. In this regime both short- and long-range 
scattering lead to the same linear density dependence of the conductivity. Our calculations imply 
that bilayer and single layer graphene have the same scattering mechanisms. We also provide an 
upper limit for the effective, density dependent spatial extension of the scatterers present in the 
experiments. 
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I. INTRODUCTION. 

Single- and bilayer graphene have attracted tremen- 
dous attention since their discover j*i due to their fascinat- 
ing and in many respects unique properties (for a review 
see, e.g., Refs. [2l-[j). Presently, one of the central issues in 
graphene research is to identify the scatterers that dom- 
inate the conductivity a. This is not only a fundamental 
question, but also the prerequisite for progress regard- 
ing the quality of the electronic properties of graphene 
devices. Since the character of the scatterers manifests 
itself directly in the electron density (n) dependence of 
(7, this function plays a key role in the corresponding 
experimental and theoretical investigationsj^i^ Most ex- 
periments in single-layer graphene (SLG) show a depen- 
dence close to (7 oc n^^, while some sub-micromcter 
extra-clean suspended samples show a square-root de- 
pendence tJ oc n^/^ consistent with the ballistic trans- 
port regime^"— . Within the standard Boltzmann trans- 
port theory in the Born approximatioufii"— the linear 
density dependence emerges from long-range Coulomb 
scattering which implies that charged impurities in the 
substrate are the dominant scatterers. Deviations from 
a linear dependence were attributed to effects of short 
range impurit ies^"— ' ^ ^ i ^ ^ , for which the Boltzmann the- 
ory predicts a = consiiii"— . This point of view has 
been challenged by theories going beyond the first Born 
approximationi^^— as well as numerical simulationsi^i^ 
which show that short-range disorder leads to a linear 
density dependency similar to the one for the long-range 
potential. At the same time, a growing evidence from 
recent experiments^i"— points to short-range impurities 
as the main scattering mechanism in SLG. 

The situation is much less clear in bilayer graphene 



(BLG). In most studies the BLG conductivity exhibits 
a linear density dependence very similar to that one of 
SLGiii2S'24 A supcrlincar dependence was reported by 
Zhu et ali^, whereas a sublinear behavior has been ob- 
served in sub-micrometer BLG samples^^ which may be 
an indication quasi-ballistic transport. In the density 
regime where these experiments have been carried out, 
the energy dispersion of BLG is parabolic. For the case 
of short-range disorder, both the standard Boltzmann 
■, as well as theories not relying on the Born 



Recent numerical 



approach* 

approximatio n^'^1^^ predict a oc 
modeling of the BLG conductivity based on the Kubo 
approach shows that sufficiently far away from the neu- 
trality point, a linear or sublinear dependence of cr on n 
is found, depending on the strength and type of of the 
disorderi^i^ As far as the long-range disorder is con- 
cerned, the Boltzmann model in the Born approximation 
predicts cr oc for unscreened, cr oc n for ovcrscreened 
and cr oc n''' with 1 < 7 < 2 for screened Coulomb 
disorde»2I. These predictions are qualitatively differ- 
ent from the corresponding predictions for SLGi^iii"— 
Therefore, in order to explain the observed density de- 
pendence (which is very similar for SLG and BLG), Das 
Sarma et al^ concluded that scattering mechanisms are 
fundamentally different for SGL and BLG, being dom- 
inated by the Coulomb impurities for the former and 
short-range impurities for the later. This is a rather sur- 
prising conclusion because scattering mechanisms in both 
systems are expected to be the same as both SLG and 
BLG are produced by the same manufacturing technique 
with the same substrate used in the measurements. 

Since conclusions regarding the nature of the scatter- 
ing mechanism in BLG are often based on the predictions 
of the standard Boltzmann approach, it is of utmost im- 
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portance to establish its range of validity and to study 
how well it models the exact conductivity. In the present 
study, we utilize the well established model of Gaussian 
disorder j^ii^iiSi^"— where the effective screening length 
^ can be chosen to describe both short- and long-range 
scattering, and derive analytical expressions for the con- 
ductivity of the BLG within the standard Boltzmann ap- 
proach in the Born approximation. We then perform 
exact numerical tight-binding (TB) Landauer-type cal- 
culations and compare them with the analytical results. 
We demonstrate that for the density regime correspond- 
ing to the parabolic dispersion the exact TB calculations 
are consistent with the Boltzmann predictions. We show 
that in the experimentally relevant regime the latter pre- 
dicts a linear density dependence of the conductivity for 
both short- and long range scattering. This density de- 
pendence is in agreement with most experiments reported 
to date. Our calculations imply that the same scattering 
mechanisms dominate in BLG and SLG. We also provide 
an estimate for the range of effective screening lengths 
for the scatterers present in the experiments. 



potential height. The correlator of this potential has also 
a Gaussian shape^^"— 

CW.i:"(r.)n..-.)^«oxp(-il), 

(3) 

where the dimensionless impurity strength is described 
by the parameter K « 40.5nimp((5/t)^(^/-\/3a)^ with 
Tiimp being the relative concentration of impurities. 

For the exact numerical calculations, we consider a 
rectangular BLG stripe of length L and width W exposed 
to the impurity potential according to Eq. ^ and attach 
it to semi-infinite leads. The conductivity of the stripe 
a{n) = ■^G(?i) as a function of the electron density n 
is obtained from the transmission coefficient T which is 
related to the conductance G via the Landauer formula 
G = T and n are computed with the aid of the re- 

cursive Green's function technique^i^. Because of com- 
putational limitations we study stripes with L/W > 1. 
The obtained results are insensitive to L/W as long as 
L/W > 1. 



II. DESCRIPTION OF THE MODEL. 

We calculate a{n) of BLG using the standard p-orbital 
nearest-neighbor tight-binding Hamiltonian 

-t^iaj.bi^i+A + h.c.) - + h.c), 

with Vi being the external potential at site i, a^^ i^tj) 
being the creation operator at sublattice A (B) and site 
in layer £ = 1,2, and the coupling integrals t = 3.16 
eV, t± — 0.39 eVi^"— In the second term the index i + A 
corresponds to the nearest neighbors to the site i. In the 
parabolic dispersion regime we follow Yuan et al^ and 
set t± = 0.5t. This effectively extends the parabolic 
band to higher densities. Without this scaling and for 
the impurity concentrations used in our numerical calcu- 
lations, the parabolic band would otherwise not be ac- 
cessible since it is essentially governed by the impurity 
band. Scattering centers are modeled by the Gaussian 
potential^^ia^-^ 

F(rO^ EV.exp(-^i^). (2) 

We refer to the case £, ~ a (appropriate for the absorbed 
neutral impurities) as short-range scattering, and to the 
case ^ ^ a (appropriate for the remote charged impuri- 
ties) as long-range scattering, with a being the C-C dis- 
tance. The potential heights are uniformly distributed 
in the range Ui S [—6, S] where 6 denotes the maximum 



III. THE BOLTZMANN APPROACH 

Within the Boltzmann approach, the conductivity 
reads^ 



a=^eVI?(^;^)(^,|?«)^ (4) 

where D{Ep) is the density of states (DOS) at the Fermi 
energy, Vp'^ is the Fermi velocity of the bilayer graphene, 
and T is the scattering time. An analytical expression 
for the Boltzmann conductivity for the bilayer graphene 
can be obtained for the limiting cases of the parabolic 
and linear dispersions, hvp\k\ <^ ^ and hvF\k\ ^ ^ 
respectively (or, equivalently, n <^ no and n ^ no, where 
no — -^fjjaj^ (^) is the critical density separating the 
parabolic and linear bands, and vp = ^ is the Fermi 
velocity of single- layer graphene). 

In this section we present analytical expressions for 
the conductivity a for the bilayer graphene. The details 
of the derivation are given in the Appendix. For the 
parabolic dispersion (n ^ no) the conductivity reads, 



m*2 K(nvp)^ (i - 1) /i(z) + (i + 1) hiz) ' ' 

( n, z <C 1 (quantum scattering); 
[ n^/^, z ^ 1 (classical scattering). 

where is the modified Bessel function, and z = 7rn^^ ~ 

(^f^) ■^ith A being the Fermi wavelength. It follows 
from the definition of z that the condition z <ti 1 corre- 
sponds to the case of quantum scattering when the Fermi 
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parabolic band 
(n < no) 


linear band 

(n > no) 


quantum scattering 
{z « 1) 


Regime I: 
G (X n 


Regime II: 
cr = const 


classical scattering 
[z » 1) 


Regime III: 
(7 oc n^/^ 


Regime IV: 

a oc n^/^ 



TABLE I. Four different regimes for the density dependence 
for the conductivity of BLG as predicted by the Boltzmann 
theory, Eqns. (O,®. It is noteworthy that the exponents for 
the linear band of BLG (regimes II and IV) are the same as 
those for SLG^^i^. 



wavelength is larger than the effective width of the po- 
tential barrier, A 3> ^, while the opposite condition z 3> 1 
corresponds to the case of classical scattering, A <C ^. 

For the linear dispersion {n 3> tiq) the conductivity 
reads. 



16e2 



h K{h{z) + h{z)) 
z <^ 1 (quanti 

z ^ 1 (classical scattering). 



(6) 



const, z ^ 1 (quantum scattering); 

^ ^ „3/2_ 



Hence, it follows from Eqs. ^ and (jH]) that the 
Boltzmann approach for the Gaussian potential predicts 
four different regimes where the density dependence of 
the conductivity a = is qualitatively different for 

parabolic and linear bands {n <C no respectively n ^ ng) 
and for quantum and classical scattering ( z ^ 1 respec- 
tively z 3> 1). The corresponding asymptotes for the 
density dependence for these four regimes are summa- 
rized in Table L 



IV. EXACT RESULTS AND DISCUSSION 

As shown in the preceding Section, for the case of the 
Gaussian potential the Boltzmann approach predicts four 
different regimes where the density dependence of the 
conductivity a = cr(n) is qualitatively different, see Eqns. 
dS]), ([U and Table L In the present Section, we compare 
these predictions with the exact TB calculations. 

Figure [1] displays the results for Regimes I and II 
(quantum scattering, z 1, for the parabolic and lin- 
ear band, respectively). Let us start with the most in- 
teresting Regime I of the parabolic band appropriate for 
most experiments (left part of the figure, n <^ ng). The 
TB conductivity exhibits a density dependency close to 
the expected linear one, cr oc n, and the TB and Boltz- 
mann conductivities are in a reasonable quantitative 
agreement. Note that for smaller scattering strengths 
{S < 2t) the density dependence of a'^^ is rather sublin- 
ear, which can be attributed to the quasi-ballistic trans- 
port regime (a transition from the ballistic to diffusive 
transport regimes in the SLG is discussed in Refs.— i^). 
The calculations presented in Fig. [T] correspond to the 
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FIG. 1. (color online) The density dependence of the TB and 
Boltzmann conductivities for Regime I (quantum scattering 
z <C 1, parabolic band n ^ no) and Regime II (quantum 
scattering z <^ 1, hnear band n 2> no); ^ = a; fx = 0.5t, 
corresponding to no = 4.4 x lO^^cm"^, for selected scatter- 
ing strengths. The TB calculations are averaged over 1000 
impurity configurations. 



short-range potential with ^ = a. Similar results are 
also obtained for a potential with a longer range, ^ = 3a 
(not shown). Because of computational limitations we 
are not able to explore BLG conductivity in the Regime 
I for larger ^. However we expect that for a given poten- 
tial strength the agreement between the Boltzmann and 
the TB conductivities improves when the potential be- 
comes smoother, i.e. when ^ increases. A similar linear 
dependence of the BLG conductivity for the Gaussian 
potential in the parameter range a ^ C ^ 5a was also 
obtained by Yuan et al^. 

Regime II of the quantum scattering case, z <^ 1, in 
the linear band corresponds to the right part of the fig- 
ure, n 3> uq. In this case the TB and Boltzmann con- 
ductivities disagree both qualitatively and quantitatively: 
the TB calculations show a linear dependence a'^^ oc n, 
whereas Boltzmann approach predicts a^°^^^ = const. 
For the case of the linear band the BLG can be effec- 
tively regarded as two weakly interacting SLG sheets. 
Therefore the SLG conductivity and the conductivity of 
the BLG in the linear band are expected to show similar 
features. Indeed, the linear density dependence of the 
TB conductivity as well as a qualitative disagreement 
between the TB and Boltzmann approaches within the 
Born approximation for z <C 1 have also been found for 
Gaussian scatterers for SLG^^. 

Why does the Born approximation for z <C 1 work 
reasonably well (in the parameter region explored here) 
for the BLG in the parabolic band (Regime I), but fails 
for BLG in the linear band (Regime II) as well as for the 
SLG? The answer to this question can be obtained from 
an analysis of the validity of the Born approximation 
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based on non-pcrturbativc approaches^SiH. Using the T- 
matrix technique for the lattice model or the partial- wave 
expansion for the continuum description one can obtain 
an exact expression for the scattering rate for the case 
of the short-range ((5-function) impurities. The condition 
for the validity of the Born approximation for short-range 
scattering was discussed by Ferreira et ali^^ and reads 



{Veff/t) kR\n{kR) < 1, SLG; (7a) 
< 1, BLG (parabolic band) (7b) 



m/tj_ 



--a, Ai 



3V3„2 
4 " 



where k is the Fermi wavelength, R = 
is the area per one C atom, c is a constant of the or- 
der of 1, and Veff is the effective impurity potential ob- 
tained using the T-matrix approach and the ab initio 
calculations^Si^. For the SLG (and thus for the BLG 
in the linear band), the condition (|7ap implies that the 
Born approximation is valid only for weak potentials, 
Veff/t ^ 5 . In contrast, condition (j7b[) is more relaxed, 
such that the Born approximation for the parabolic band 
of BLG is expected to work at least qualitatively even for 
strong reahstic potentials, Veff/t ^ 60, appropriate e.g. 
for adsorbed hydrogen^S. In the present study, we present 
results for the short range potential with Veff/t < 4. For 
larger Vef f /t the structure at hand for a given impurity 
concentrations and a system size enters the localization 
regime, where the Boltzmann approach is not applicable. 
Although we expect that the Born approximation would 
be qualitatively correct even for such strong potentials, it 
would be interesting to explore this case using numerical 
approaches capable to treat larger size structures, such 
as e.g. the Kubo approach^. 

Let us now turn to the case of classical scattering, 
z 3> 1, corresponding to Regimes III and IV (parabolic 
and linear bands respectively). We start by noting 
that for z ^ 1 the conductivity of SLG is well de- 
scribed by the Boltzmann approach within the Born 
approximationiSi^. Wc therefore expect a similar agree- 
ment between the Boltzmann and the TB calculations 
for BLG for the case of classical scattering in the linear 
band regime (i.e. in Regime IV). The conductivity of 
BLG in Regime IV (exemplified by choosing ~ 16a) is 
shown in Fig. [21 According to the expectations, the TB 
calculations are in excellent qualitative agreement with 
the Boltzmann predictions showing the expected density 
dependence a cx n'^/^ and exhibiting the decrease of con- 
ductivity as the impurity strength S is increased. How- 
ever, the magnitude of the TB conductivity differs from 
its Boltzmann counterpart by a factor 3-6 with increasing 
deviation as d increases. A similar quantitative discrep- 
ancy between the TB and the Boltzmann conductivities 
is also found for Regimes I and III. One of the reasons 
for this discrepancy is that due to computational limita- 
tions, the dimension of the structure L is not sufficiently 
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FIG. 2. (color online) The density dependence of the TB 
and Boltzmann conductivities for the Regime IV (classical 
scattering 2 3> 1, linear band n 3> no), ^ = 16a ; t± = 
0.39 eV, corresponding to no — 2.67 x 10^'^cm^^, for various 
scattering strengths. The TB calculations are averaged over 
1000 impurity configurations. The sample size is L x W = 
320 X 20nm^ 



large to achieve a truly classical diffusive regime with 



L > Ztr > h 



(8) 



where L 



VT 



is the mean free path and Ip 



is the Fermi wavelength. For the scattering strengths 6 
used in calculations typical values of the mean free path 
are kr ~ L/2 - L/10 (L = 320 nm in Fig. [2]and L = 160 
nm for Figs Ill3p . and ~ 10 nm. We expect that the 
agreement between the TB and Boltzmann conductivi- 
ties would improve for larger structures where the con- 
dition Eq. ([5]) is better satisfied. Another reason for the 
discrepancy can be related to the utilization of the Born 
approximation, the applicability of which improves as the 
electron energy increases. 

We finally turn to Regime III (classical scattering, 
2; ^ 1, in the parabolic band). It is exemplified in Fig. 
[3] by choosing ^ = 32a. The TB calculations are in a 
qualitative agreement with the corresponding Boltzmann 
predictions showing the density dependence a cx n^/^ ex- 
pected for this regime and exhibiting a decrease of con- 
ductivity as the impurity strength S is increased. We con- 
clude this discussion by summarizing that for the classical 
scattering, z 1, the TB calculations are consistent with 
the Boltzmann approach for both parabolic and linear 
bands, whereas for the quantum scattering, z <C 1, the 
TB and Boltzmann calculations agree for the parabolic 
band and disagree for the linear band. 

Let us relate our findings to available experimental 
data. Practically all the experiments on BLG struc- 
tures reported so far have been performed for relatively 
low electron densities where the dispersion relation is 
parabolic. In most of these studies the BLG conductivity 
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FIG. 3. (color online) The density dependence of the TB and 
Boltzmann conductivities for the Regime III (parabolic band, 
n <C no and weak scattering z ^ V). ^ — 32a; t_i_ = 0.5t 
(corresponding to no = 4.4 x 10^'' cm^^), shown for different 
scattering strengths. The TB calculations are averaged over 
1000 impurity configurations. 

exhibits a linear density dependence i^i^i^i. Out of four 
regimes considered above, this is consistent only with the 
Regime I (quantum scattering in the parabolic band) pre- 
dicting a linear density dependence for both short- and 
long range scatterers. We therefore conclude that Regime 
I is the regime corresponding to realistic BLG samples, 
and the experimentally observed linear dependence of the 
BLG conductivity can be caused by both short- and long- 
range scattering. Note that exact TB numerical calcula- 
tions for SLG for the case of quantum scattering, z <C 1, 
also give the same linear density dependence for both 
short- and long-range Gaussian scatterersi^. As the cal- 
culated and experimental conductivities of SLG and BLG 
are essentially the same, this strongly suggests that the 
scattering mechanisms for SLG and BLG can not be dif- 
ferent. 

As far as the classical scattering z 1 in the parabolic 
band is concerned (Regime II), the Boltzmann theory 
and the TB numerical calculations predict a supcrlin- 
ear density dependence a a n^^^. This dependence has 
never been observed in BLG structuresi^ This leads us 
to conclude that the regime of the classical scattering, 
z ^ 1, is not accessible experimentally for BLG. (It is 
noteworthy that a similar conclusion, namely that the 
regime of classical scattering z 3> 1 is not achieved in 
SLG, was reported in Ref.[13). Therefore, the condi- 
tion z = 7rn^^ ~ 1 provides an upper limit for the ef- 
fective, density dependent spatial extension of the scat- 
terers present in the experiments. For, example, for a 
typical electron density n = lO^^^cm"^, we obtain ^ ~ 5.6 
nm (~ 22 lattice constants). This value of the effec- 
tive screening length ^ is consistent with the result re- 



ported by Ghaznavi et al^ for the Poisson-Thomas- 
Fermi screening of a charged impurity where a potential 
drops by a factor of 10 at distances 1-5 nm. 



V. CONCLUSIONS 

We have studied the conductivity of BLG utilizing a 
model of Gaussian disorder where the effective screen- 
ing length ^ can be chosen to describe both short- and 
long-range scattering. Analytical expressions for the con- 
ductivity of the BLG within the standard Boltzmann ap- 
proach in the Born approximation have been derived and 
compared to exact numerical tight-binding Landaucr- 
type calculations. Our main findings can be summarized 
as follows. 

• The Boltzmann approach for the Gaussian poten- 
tial predicts four different regimes where the den- 
sity dependence of the conductivity a = <j{n) 
is qualitatively different for parabolic and linear 
bands {n <^ uq respectively n 3> n-o) and for the 
quantum and classical scattering ( z <C 1 respec- 
tively z ^ 1), see Table I for the summary. For the 
classical scattering, z ^ 1, the TB calculations are 
consistent with the Boltzmann approach for both 
parabolic and linear bands, whereas for the quan- 
tum scattering, z <C 1, the TB and Boltzmann cal- 
culations agree for the parabolic band and disagree 
for the linear band. We discuss and analyze the ap- 
plicability of the Born approximation for BLG and 
compare it with the case of SLG. 

• By comparing the obtained density dependencies 
of the conductivity with available experimental re- 
sults we conclude that realistic BLG samples are 
in the regime of the parabolic band, n uq, and 
of quantum scattering, z <C 1. Here, both short- 
and long-range scattering lead to the same linear 
density dependence of the conductivity. We also 
compare the conductivities of SLG and BLG and 
conclude that they have the same scattering mech- 
anisms. 

• An upper limit for the effective, density dependent 
spatial extension of the scatterers present in the 
experiments is provided. For, example, for a typical 
electron density n = lO^^^cm"^, we obtain ^ — 5.6 
nm (corresponding to 22 lattice constants). 
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Appendix A: Calculation of the Boltzmann 
conductivity a for bilayer graphene 



In this appendix we present a derivation of the Boltz- 
mann conductivity for bilayer graphene for the case of the 
Gaussian potential. According to Eq. , calculations of 
the conductivity a requires knowledge of the scattering 
time T, the DOS at the Fermi energy D{Ep), and the 
Fermi velocity Vp*^ of the bilayer graphene. 



1. Basics of bilayer graphene 

In this section we briefly summarize basic properties of 
bilayer graphene including the dispersion relations, wave 
functions, DOS and the Fermi velocity that will be sub- 
sequently used below to calculate the conductivity a and 
the scattering time t. 

We write the four component wave function of the bi- 
layer graphene in the form i/-" = (c"^S , c^^ , c^^ ) e**^"" 
where A and B correspond to the A and B sublattices, 
and 1 and 2 to layers, and Q is the wave vector— i^. Substi- 
tuting this wave function into the Schrodinger equation 
with the Hamiltonian Eq. ([T]) and expanding Q in the 
vicinity of the K point, Q = K + k. wc obtain 



( ' 


./(q) 




\ 








/*(q) 













= -E 












r(q) 






V 





/(q) 


/ 









where E = E/t, = t±/t, and 



/k = {kx 



... 3a , hvpk 



) = -|./k|e-^ 
9, ^tan^. 



(Al) 
(A2) 



The eigenvalues of Eq. (jAip are easily obtained. 



Sl 32- 




, si,s2 = ±1, (A3) 



with the corresponding eigenfunctions 

/ -A \ 

SlS2-Bk 
\-SiS2,fkJ 



Mr) = 



(A4) 



where the normalization constant C = 2(i?^ + |/kP). 

The dispersion relation Eq. (|A3() is illustrated in Fig. H) 
It is instrumental to consider two limiting cases where 
the dispersion relation Eq. (|A3|) shows qualitatively dif- 
ferent behavior, (a) low energies, |/k| < {hvF\k\ < 




si = -l 



\ / 

hvp\k\ = V2tj_ 



FIG. 4. A schematic illustration of the dispersion relation for 
a bilayer graphene, Eq. (|A3|) . 



iij^) and (b) high energies, |/k| < {hvF\k\ > ^tj_). 
Using the relation between the wave number and the elec- 
tron density in two dimensional (2D) electron systems, 
k = y/im, the conditions (a) and (b) can be rewritten 
in terms of the electron density as follows, (a) low den- 
sities, n ^ no, and (b) high densities, n ^ ng, where 



no 



^(^^^yi it±/t) is the critical density separating the 



parabolic and linear bands. We will consider only the 
conduction band (si = 1); the same expressions for 
the conductivity can be obtained for the valence band 
(si = -l). 



a. low densities, n ^ no (parabolic band) 



Introducing the effective mass m* = -^r^ 



, the disper- 
3]) for the lowest band {s2 = —1) can 



sion relation Eq 
be re- written in a form familiar for conventional semicon- 
ductor systems, 



E = 



2m*' 



(A5) 



Hence, the Fermi velocity for the bilayer graphene Vp'^ 
and its DOS have the same form as the ones for conven- 
tional 2D electron systems. 



,BG 



= hk/m*- D^^(E) = 



2m* 



(A6) 



where the factors gy = gs — 2 account for the valley and 
the spin degeneracy. 



b. high densities, n S> no (linear bands) 

For the high densities the dispersion relation consists 
of two linear bands (as in the case of the the single layer 
graphene) separated by the energy interval t±, 



7 



E ; 



hvF\k\, S2 = —1; 

/lWF|fc| + t±, S2 = 1. 



(A7) 



Each of these bands has the same dispersion as the one 
of single- layer graphene (SLG); therefore, in this case 
Vp*^ = I'F, and the DOS reads 



1SG^E^^ _ onSLG/pn^i ^ 9v9sE ^ AE 



D^^{E) = 2D^^^{E) 



(A8) 



where a factor of two accounts for two bands. 



2. Calculation of the scattering rate 



Within the Boltzmann approximation the scattering 
rate reads^i 



where in the Born approximation Wk.k' is given by 
Fermi's Golden Rule, 

W^k.k' = ^ \V^,M'\'KEk - Ek'), (AlO) 
n 

with 9 being the angle between the wave vectors k and k' , 
and Vk k' is the matrix element of the scattering potential 



'4^iVir)ij.,dr. 



(All) 



Using the wave function Eq. (jA4p we obtain for the ma- 
trix element, 



2(EkEk> + fkrk. + rjk' 



2dEl + \h\^jEl + \h 



:f/„ (A12) 



where Uq = J V{r)e'''^'^dr is the Fourier transform of 
the scattering potential, and q ~ |k — k'| = 2fcsin|. 
We calculate |C/qP making use of the Wiener-Khinchin 
theorem^ stating that the square modulus of the Fourier 
transform of a function is equal to the Fourier transform 
of its correlator, |J7gP = / e"i''C(r)dr. Using the expres- 
sion Eq. ^ for the correlator of the Gaussian potential, 
we get 



(A13) 



One can obtain analytical expressions for the scattering 
rate Eq. (|A9p for two limiting cases of the low densities 
{n ^ no) and the high densities (n ^ no) where elec- 
trons obey respectively quadratic and linear dispersion 
relations. 

a. low densities, n <^no (parabolic hand) 
Using Eqns. (fJ2|l and (|J5|) in Eq. (|A12|) we obtain 



1 + cos 261 



U, 



(A14) 



1 f dk! 

-=y (^^k,k'(l-cos0), (A9) i_ ^D[E) 



The expression for the scattering rate Eq. (jA9p reads, 

- = / d0(l + cos20)(l-cos0)|[/j2, (A15) 

T 4/i Jo 

where we used the expression Eq. (jA6|) for the DOS. The 
Fourier transform of the Gaussian potential is given by 
the expression Eq. (|JT3)) . Substituting it into Eq. (|JT5|) 
and performing integration utilizing that q = 2k sin |, we 
obtain for the scattering rate 



(A16)" 

where z = fc^^^ = 7rn^^, and Ii,{z) are the modified 
Bessel functions of the order ly. Substituting this expres- 
sion into Eq. Q and using the definition of the DOS and 
Vp'-^ for the bilayer graphene, Eq. (jA6p . we finally obtain 
Eq. ([S]) for the Boltzmann conductivity. 



6. high densities, n ^ no (linear bands) 
Using Eqns. ^M^l and (Mil in Eq. (|Il2|) we obtain. 



1^4,1 



\Ua 



(A17) 



Substituting the expression for into Eq. (jA9|) and using 
Eq. (|A8p . we obtain for the scattering rate. 



1 D 



BG 



T I6h Jq 



d9 (l + cos 9 f {1- cos e)\Ugf. (A18) 



Substituting Eq. (|A13[) into this equation and perform- 
ing integration using that q = 2fcsin |, we obtain for the 
scattering rate. 



1 ttD 



BG 



.K{hvF?e-^ [ IM±1M. ) . (A19) 



Substituting this expression into Ec^. (|4]) and using the 
definition of the DOS and Vp'~^ for the bilayer graphene, 
Eq. (|A8|) . we finally obtain Eq. ^ for the Boltzmann 
conductivity. 
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